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Abstract 

In this work we define a formal notion of a quantum phase crossover for certain 
Bethe ansatz solvable models. The approach we adopt exploits an exact mapping 
of the spectrum of a many-body integrable system, which admits an exact Bethe 
ansatz solution, into the quasi-exactly solvable spectrum of a one-body Schrodinger 
operator. Bifurcations of the minima for the potential of the Schrodinger opera- 
tor determine the crossover couplings. By considering the behaviour of particu- 
lar ground-state correlation functions, these may be identified as quantum phase 
crossovers in the many-body integrable system with finite particle number. In this 
approach the existence of the quantum phase crossover is not dependent on the 
existence of a thermodynamic limit, rendering applications to finite systems fea- 
sible. We study two examples of bosonic Hamiltonians which admit second-order 
crossovers. 

Keywords: quantum integr ability (Bethe ansatz), Bose-Einstein condensation (the- 
ory) 



1 Introduction 



Quantum phase transitions may occur in the ground state (i.e. at zero temperature) of 
quantum systems as an interaction coupling is varied. Such phase transitions may be 
thought of as being driven by quantum fluctuations, in analogy with thermal fluctuations 
underpinning thermal phase transitions. In fact in many cases there is a correspondence 
between a quantum phase transition in d dimensions and a thermal phase transition in 
(d + 1) dimensions [I]. One of the important mathematical tools in the study of thermal 
phase transitions is the thermodynamic limit, where the number of particles is taken to 
infinity. In this limit phase transitions can be associated with discontinuities in certain 
physical quantities derived from the free energy. 

For finite (i.e. mesoscopic) systems we cannot appeal to this notion of discontinuities 
to characterise a sharp change between quantum phases, as the transitions are smooth. 
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Nonetheless for many finite systems there is a sense of a crossover between different 
quantum phases, and it useful to be able to characterise such a crossover, beyond saying 
that it is what occurs in a finite system in cases where there is phase transition in the 
thermodynamic limit. Quantum phases in finite systems have previously been studied 
in [2] in the context of the Interacting Boson Model (IBM). There it was argued that 
quantum phase transitions can be identified via "shape transitions" of an effective po- 
tential energy surface, defined in terms of classical variables, which becomes exact in the 
thermodynamic limit. Our goal here is to formulate an approach which entirely avoids 
the use of the thermodynamic limit, for applications to cases where taking this limit may 
not be desired or justified. For example there are known models where crossovers be- 
tween ground-state phases only exist in the finite case. One possible situation is that the 
boundaries in parameter space between such ground-state phases may merge together as 
the thermodynamic limit is approached. An example of this is seen in the classical field- 
theoretic analysis of the attractive non-linear Schrodinger equation, where the transition 
coupling between the uniform regime and the broken symmetry soliton regime occurs at 
a coupling which scales as iV _1 for large N, where N is the number of particles [5]. In the 
thermodynamic limit the transition coupling degenerates to zero (free- field case). For the 
Dicke model there are several "critical points" which degenerate in the thermodynamic 
limit [4J. Another problematic scenario arises if taking a particular limit of a coupling 
parameter does not commute with taking the thermodynamic limit, as is known to hap- 
pen in the weak coupling regimes of the BCS model E] and the Bose gas with delta 
function interactions [7J. Moreover for some systems, such as the attractive case of the 
Bose gas, there are added technical difficulties in defining the thermodynamic limit, as 
the ground-state energy per particle is not finite in this limit [7J [8]. 

Motivated by the above considerations, below we give a characterisation for a quantum 
phase crossover (QPC) in a one-body system with an external potential, for which there 
is no notion of a thermodynamic limit. Our definition for a QPC is given in terms of 
properties of the corresponding classical system. We then illustrate how this result can 
be used to investigate QPCs in Bethe ansatz solvable many-body interacting systems 
with finite particle number, through a manner which avoids taking the thermodynamic 
limit. The key to this approach, as will be detailed below, is to exploit the Bethe ansatz 
solution to perform an exact mapping from the spectrum of the many-body interacting 
system into the spectrum of a one-body system in a potential. 

2 Quantum phase crossovers 
2.1 One-body systems 

We start with the Schrodinger operator (SO) eigenvalue equation in one dimension 



where for the potential V(x) it is assumed there is a bifurcation of the global minimum 
at Xq, such that we may consider the approximation 



with Vi > 0. A classical treatment of the above is equivalent to Landau theory, as used 
in the study of thermodynamic phase transitions. This approach has been previously 



+ V(x)ip k = E k ip k 




dx 2 



V(x) &V - 2V 1 (x - x ) 2 + V 2 (x - x ) 4 
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discussed in [2J. The classical ground-state energy E is given, to leading order in Vi, by 

P J V forVi<0 

~ | Fo _ V 2 /V2 for Vi > Q (2) 

Using this result as an approximation to the ground-state energy for the quantum sys- 
tem, we can appeal to the Hellmann-Feynman theorem to approximate the ground-state 
position fluctuations: 

/( . 2X _ ldE f forVi<0 

\{X — Xq) ) - 



2 dV l 1 Vx/V 2 for > 0. 

In a full quantum analysis we must expect that quantum fluctuations will smooth out 
the discontinuity in the derivative of ((x — x ) 2 ) (and in particular ((x — x ) 2 ) > for all 
Vi), so there is no quantum phase transition in the traditional sense. Nonetheless, the 
above classical analysis indicates that around the crossover coupling V\ = we should 
expect a sharp change in the behaviour of ((x — x ) 2 ) indicative of a crossover between 
a localised state and a Schrodinger cat state. Therefore we formally define a QPC for a 
one-dimensional SO as follows: 

Definition 1 Consider the SO equation where the potential V(x) smoothly depends 
on a dimensionless coupling parameter 7 . Treating as a classical problem, approximate 
the ground-state energy as the minimum of the potential via Eq = mmV(x). If m is the 

smallest integer for which d m Eo/d , y m is discontinuous at some coupling 7 C; we say there 
is an mth-order QPC of the quantum system at j c . 

Of interest to the examples considered below are second-order crossovers, for which 
d 2 E /d'j 2 is discontinuous. 



2.2 Finite many-body systems 

In the thermodynamic limit there are known examples for which there is a one-to-one 
correspondence, established via the Bethe ansatz, between the spectrum of an integrable 
model (IM) and the spectrum of a SO 0, QUI [EE]. For IMs acting on finite-dimensional 
Hilbert spaces an analogous scenario can hold, where the spectrum of the IM maps into 
the spectrum of a SO, the so-called quasi-exactly solvable (QES) sector [T2],[T3]. Such a 
mapping is obviously not one-to-one, but typically the QES spectrum of the SO corre- 
sponds to low lying energy levels. In certain cases, such as the examples studied below, 
the mapping is faithful between the ground-state energies of the IM and the SO. This fact 
provides the means to explore QPCs in the IM. 

Let H = Hi^f) denote the Hamiltonian for some IM acting on a finite-dimensional 
Hilbert space where 7 is a dimensionless coupling parameter. Let \^fj) denote the eigen- 
states of H with energy levels Sj and the convention that Sq is the ground-state energy. 
Now assume that the energies Sj can be mapped to energies E^ of a SO such that for each 
j there exists a k satisfying 

£ 3 = xE k (3) 

for some positive scale factor x which is independent of 7, j and k. Since the mapping 
of the spectrum of the IM to the corresponding SO is not one-to-one, it is important to 
assert that the mapping is faithful between the ground-state energies. In our examples, 
following the arguments in [T^J [T3], this will be guarranteed by the following oscillation 
theorem |14j : 
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Theorem 1 Consider a SO with locally bounded potential V(x) satisfying V(x) —>■ oo as 
\x\ — > oo. Let ipk, k = 0, 1, . . . oo denote the eigenf unctions of the SO with eigenvalues E^ 
respectively, ordered such that Ej < Ef. whenever j < k. Then ipk has precisely k (real) 
zeroes. 

A corollary of the theorem is that the ground state wavefunction of the SO has no zeroes. 

Now assume that the mapping ([3]) does hold between the ground state energies. Defin- 
ing the operator A = dH/d^ which acts on the Hilbert space of the IM, we again use the 
Hellmann-Feynman theorem to deduce that for the ground state 

M\ = d£o = ^ dEp _ ^ dE 
$7 (97 (97 

Thus the behaviour of the ground-state correlation function (A) will display a sharp 
change if d 2 E /dj 2 is discontinuous. This leads to 

Definition 2 Suppose the ground state energy £ of an IM exactly maps to the ground 
state energy E of a SO through ^). If the SO exhibits a second-order QPC at some 
dimensionless coupling 7 C as in Definition 1, then we say that the IM also exhibits a 
second- order QPC at j c . 

Having outlined the theroetical considerations, we now apply it to specific models. 



3 Examples 

3.1 Atomic-molecular bosonic model 

The first model we look at describes the interconversion of bosonic atomic and di-atomic 
molecular modes [15]. The Hamiltonian is 

H = -n a + ^(aVfr + feW) (4) 

where o) and W denote the creation operators for atomic and molecular modes respectively 
and as usual n a = a^a, = b'b. The total particle number N = n a + 2^ is conserved. 
In addition the Hamiltonian is invariant under the transformation (a, a') — > (—a, —a'). 
Since the change Q — > —Q is equivalent to the unitary transformation (b, &') — > (—b, —tf) 
we restrict to Q > 0. 

The Hamiltonian admits an exact Bethe ansatz solution [16J. The solution depends 
on a discrete variable p which takes values or 1, depending on whether the total particle 
number is even or odd. The exact solution gives the energy levels as 

M 

J3 = tf(M+!)+nX;«i (5) 

3=1 

where the parameters {vj} are the roots of the Bethe ansatz equations 

Above, 7 = 5/Q is defined as the dimensionless coupling, the number of roots M is related 
to the total particle number through N = 2M +p, and the dimension of the Hilbert space 
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Figure 1: Behaviour of the coherence correlator ([7]) for the ground state of (T4]). The 
curves shown are for total particle number N = 20,30,40. The inset shows (n a ) /N, 
the average fractional occupation of unbound atoms in the ground state. Below the 
crossover coupling S/(ri\^N) « 1.4 the ground state is a coherent superposition of atomic 
and molecular states, while above the crossover coupling the ground state predominantly 
consists of molecular bosons. 



is M + 1. From the ground-state energy, one can use the Hellman-Feynman theorem to 
compute the ground-state correlations 
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where 

9 = - (oVfe + tfaa) /2 (7) 

is the coherence correlator. 

Results of numerical analysis of the exact solution are shown in Fig. [TJ which has 
been taken from [16] . Despite the very small number of particles, it can be deduced that 
below the crossover coupling 8/(Tly/~N) ~ 1.4 the ground state is a coherent superposition 
of atomic and molecular states, while above the crossover coupling the ground state pre- 
dominantly consists of molecular bosons (see also Fig. 1 of [IT]). The crossover sharpens 
with increasing N but becomes singular in the thermodynamic limit iV — > oo, implying 
Q — > and the eigenstates approach Fock states. However in this limit the ground-state 
energy per particle is not finite when S/Q < (see comments in the Conclusion). We note 
that the qualitative features of (n a ) /N shown in the inset are the same as those shown 
in Fig. 3 of [2] for the IBM. There it was argued the result could be interpreted as a 
second-order crossover, even for particle numbers of the order 10 1 . For (j3J) we will show 
that this same conclusion can be obtained from our approach described above. 

To derive the value of the crossover coupling from the Bethe ansatz solution, we map 
the spectrum into that of a SO with a sextic potential. The procedure can be found in 
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Figure 2: Generic behaviour of the sextic potential ([9]) for different values of the dimen- 
sionless coupling 7, where the crossover coupling is given by 7 C = \/2N + 3. Here we 
take N = 100. For 7 > 7 C (solid line - 7 = 30) there is a single minimum at x = 0. For 
~lc < 7 < 7c (dash line - 7 = —10) or 7 < — 7 C (dot line - 7 = —30) there are two global 
minima at ±Xq given by ffTUl) . 



[T2] . so we simply give the result (see also [13]). For a given solution of (jEJ) with energy 
given by (JSJ), we set 



(2 4 \ ^ / 2 



Then ip(x) satisfies (pQ) with \ — ^ an d potential 

v ; 4 16 32 256 V ' 

The number of QES states (jSj) is the dimension of the Hilbert space for M + 1. It 
is apparent that the states (|HJ) are even (odd) functions of x, with an even (odd) number 
of zeroes for p — (p = 1) and at most 2M (2M + 1) zeroes. In view of Theorem 1 
no two states can have the same number of zeroes. For p = the QES states must be 
the states ipk such that k = 0, 2, 4, ... , 2M. Therefore the ground state of the SO with 
potential lies in the QES sector, and the mapping of the spectrum of (j4j) into that of 
the SO is faithful for the ground-state energies. For a system with odd particle number, 
the QES states correspond to the k — 1,3, ... , 2M + 1 states of the associated SO. For 
this case one can alternatively consider the potential to be restricted to the half line x > 
with the hard-wall boundary condition V(0) = 00, requiring ^(0) = 0. Then the QES 
states have an even number of zeroes in (0, 00) and correspond to the states ipk such that 
k — 0,2,4, ... , 2M, which includes the ground state. 
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For 7 > and (7 2 — 3 — 2iV) > the potential attains its minimum at Xq = 0, 
while for ( r y 2 — 3 — 2iV) < or 7 < minima occur at 

^ = ^v / 47 2 - 3( 7 2 - 3 - 2iV)-^. (10) 

This analysis identifies the crossover coupling as 7 C = \/2N + 3. As 7 — * (7 C )- we find 
xq ~ ±(\ / 2N + 3 — 7) 1 / 2 . The above shows that for the potential (EJ) there is a second- 
order QPC of the same universality class as the Landau (or mean-field) theory, which 
is also the same class for the second-order crossover of the IBM [2]. In the classical 
approximation the behaviour of the correlation functions is found to be 

- C ~ 7 - 7 C , 7 -> (7 C )_ 

e - e c ~ 0, 7 - (7c)+, 

with # c = 0. Similarly results hold for (n a ). These are consistent with Fig. [TJ 
3.2 Attractive two-site Bose— Hubbard model 

Quantum tunneling of bosons, based on a two-mode approximation, can be described by 
the two-site Bose-Hubbard Hamiltonian [18] : 

H = - k -{n 1 -n 2 f- £ -{b\b 2 + b\b l ) (11) 

where bp j = 1,2 denote the single-particle creation operators associated with two bosonic 
modes and rii = b\b\ and n 2 = b\b 2 are the corresponding number operators. The total 
particle number iV = rii + n 2 is conserved. We only consider the attractive case for which 
k > 0, and take £ > 0. The case S < can be obtained by the unitary transformation 

(h, b{) (-&!, 

To derive the Bethe ansatz solution we follow the approach of [131 EE] • We start with 
the Jordan- Schwinger realisation of the su(2) algebra generators 

S+ = b\b 2 , S~ = blh, S^^Nt-Nz) 

which satisfy the commutation relations 

[S z , S±} = ±S ± , [S+, S-] = 2S Z . (12) 

The realisation is (N + l)-dimensional when the constraint of fixed particle number iV = 
Ni + N 2 is imposed. In terms of this realisation the Hamiltonian may be written 

H = ~ k -{S z f- £ - (S+ + S-). (13) 

The same (N + l)-dimensional representation of su{2) is given by the mapping to differ- 
ential operators 

d N Q+ AT 2 d C - d 

S =u- — , S + = Nu-u— , S = — 

au 2 du du 
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acting on the (N + l)-dimensional space of polynomials with basis {l,u,u 2 , ...,u N }. We 
can then equivalently represent (|T3|) as the second-order differential operator 

ku 2 d 2 l, I/Ar 1N ^ 2U d fciV 2 SNu 

+ - - 1)« - 5(1 - M 2 )) — - — — . (14) 



2 dw 2 2 v v y v y/ dw 8 2 

Solving for the spectrum of the Hamiltonian fill I) is then equivalent to solving the eigen- 
value equation 

HQ = EQ (15) 

where H is given by ( f!4l) and Q(u) is a polynomial function of u of order iV. To obtain 
the Bethe ansatz solution, we first express Q(u) in terms of its roots {vj}: 

N 

Evaluating (fl5l) at u — v k for each k leads to the set of Bethe ansatz equations 

g(i-,a+*(i-*H = f _g_ fc= 

v 

Writing the asymptotic expansion Q(u) ~ « w — u^ 1 and by considering the terms 

3=1 

of order iV in ffTBI) . the energy eigenvalues are found to be 

3=1 



We may exactly map the spectrum of (TTTi) into that of the SO equation (jTj). Setting 
j — £/k, x = k/2 and 



iV 



■^(x) = exp(— 7 cosh(x)) Y\ ( ex P y ~z ) ~~ ^ ex P ( "o" 

then, as a result of (fT5l) . ^(x) satisfies (JTJ) with the potential 

= 7 2 sinh 2 (x) - (A^ + l) 7 cosh(x) (18) 

whenever the {vj} satisfy (fl~6l) . Using the same argument as in the previous example, it 
can be established that this mapping is faithful for the ground-state energies. 

The potential (TTBl is a double Morse potential, which has previously been studied as 
a quasi-exactly solvable potential [131 [20]. It is a single well potential when 27 > (N + 1) 
with a mimimum at Xq = 0, and a double well for 2 7 < (N + 1) with minima at Xq = 
± cosh _1 [(iV + l)/27]. We identify the crossover coupling as 7c = (N + l)/2, and deduce 
that as 7 — > (jc)- 

±x ~((iV + l)/2-7) 1 / 2 
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establishing that (jTTj) is in the same universality class as (jl]). The coherence correlator in 
this instance, 6 = — (&{&2 + b\bi) /2, displays the classical critical behaviour 

6»-6» c ~7-7 c , 7 -> (7J _ 
- # c ~ 0, 7 - (7c)+ 

where 9 C = k^ c . We note that (1181) . with change of variable x — > ix, was previously 
derived in [2T] as an exact quantum phase model for the repulsive case of (jl II) . 

The QPC coupling 7 C = (iV+l)/2 agrees, to leading order in N, with the results of [22] . 
In that work the QPC arises as the onset of broken symmetry in the mean-field analysis. 
The result is also in agreement with [23], where the QPC was identified by a numerical 
analysis of wave function overlaps (cf. [24]) and ground-state entanglement measured 
by the von Neumann entropy. In particular it was found that there was a peak in the 
ground-state entanglement at the crossover coupling. This is consistent with the claim 
of [25] that a peak in the ground-state entanglement occurs when there is a supercritical 
pitchfork bifurcation of the global minimum of the Hamiltonian in the phase space of the 
semi-classical limit. For (jTTj) , the semi-classical dynamics have been studied in [26] and it 
is seen from the equations of motion that such a bifurcation does occur at 7 C (to leading 
order). It is interesting to note that in contrast there is no peak in the entanglement at 
7 C for the Hamiltonian (j3J) [27J. For this latter case it has been shown in [28J that the 
bifurcation occuring in the semi-classical phase space is not of the supercritical pitchfork 
type. 

4 Conclusion 

We have developed a technique for determining QPCs in certain finite IMs admitting a 
Bethe ansatz solution, and applied it to the Hamiltonians (jTjfTTj) . In both cases we iden- 
tified a dimensionless crossover coupling j c , whose role in defining the QPC is supported 
by the behaviour of particular ground-state correlation functions. The significance of the 
results is reflected by the fact that in both cases 7 C is iV-dependent, which highlights the 
potential difficulties in studying quantum phases in the thermodynamic limit. Further- 
more, for the Hamiltonian (jlj) the ground state energy per particle Sq/N scales as N l l 2 
whenever 7 < (as can be deduced from equations ([H]) and (TTUT) ) . For the Hamiltonian 
(jTTj) the ground-state energy per particle at 7 C is given, in the classical approximation, by 
Sq/N = kN(l + l/N) 2 /4. In both cases the ground-state energy per particle is not finite 
in the thermodynamic limit. 

Both models we have studied have only two degrees of freedom. In generic bosonic 
systems it is quite reasonable for a system to have a low number of degrees of freedom, 
since the bosonic statistics allow for arbitrarily high particle number independent of the 
degrees of freedom. This situation is in stark contrast to fermionic systems where, due 
to the exclusion principle, increasing the number of particles is associated with increasing 
the number of degrees of freedom when taking the thermodynamic limit. A consequence 
of this difference between bosonic and fermionic systems is that for bosonic systems with 
small number of degrees of freedom, such as the two studied here, it is not possible to apply 
renormalisation group techniques to determine the nature of the ground-state phases. 

Finally, we comment on the feasability of extensions to other models. Mappings from 
the spectrum of a many-body Hamiltonian to an SO can be constructed in cases where 
there is a spectrum generating Lie algebra such that the Hamiltonian commutes with the 
Casimir invariants. The su{2) case has been discussed in [2U], where it was argued that 
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31 classes of quasi-exactly SOs exist. Other examples can be found in [13J. However it 
remains a challenge to determine the full extent of applicability of the methods we have 
introduced here. 
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